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1. Introduction 

In the last years there has been a growing interest in nonequilibrium phase transitions 
[H |2]. The absence of a general theory nonequilibrium systems originates to a number 
of open problems, even in one-dimensional systems that, in the equilibrium regime, 
generally are exactly solvable. Usually, numerical simulations are a userful technique 
in the study of phase transitions and critical phenomena, and its power has been 
growing with the increasing capacity of the computers and the development of new 
simulational techniques. Non-equilibrium models are particularly suited for simulations. 
However, other approximate approaches may be complementary in the study of these 
phenomena. Thus, it is of interest to study the systems by other techniques. One of 
the successful techniques are the series expansions |3J, which in some cases lead to very 
precise estimates for the critical properties that characterize these transitions [31 Hj. 

Among the non-equilibrium models, are the ones that exhibit absorbing states, 
that is, states that may be reached in their dynamics, but transitions leaving them are 
forbidden. Since such models do not obey detailed balance, they are intrinsically out of 
equilibrium. The most studied system for this class of problems is the so called contact 
process (CP) model [5], a toy model for the spreading of an epidemic. This model 
displays a transition between an absorbing and an active state with critical exponents 
belonging to the directed percolation (DP) universality class (6]. In addition, the CP 
model is related to the Schlogl's lattice model for autocatalytic reactions [7] and the 
Reggeon Field Theory [8]. 

Many variants of this model have been studied [9j [101 [TTJ, [TJ], most of them 
belonging to DP class also. In fact, the robustness of this universality class is an evident 
characteristic of these models. Such robustness is explained by the conjecture that all 
models with phase transitions between active and absorbing states with a scalar order 
parameter, short range interactions and no conservation laws belong to this class [13] . 

One of these variants is the CP with diffusion [14] , which exhibits a critical line. 
This line begins at the critical point of the model without diffusion and ends in the 
infinite diffusion rate limit, where the critical properties of the system approaches those 
predicted by the mean-field approximation. The mean-field behavior of the model in 
the limit of infinite diffusion rate may be understood considering that, since diffusion 
processes are dominating the evolution of the system in this limit, creation processed 
are effectively determined by the mean densities, as is supposed in the mean-field 
approximation. This change of behavior at the infinite diffusion rate limit, between the 
critical behavior of the DP class and the one predicted by the mean-field approximation, 
characterizes a crossover of the critical properties in the neighborhood of a multicritical 
point. As in the equilibrium case we may then write any density variable, in the 
neighborhood of a multicritical point, as the following scaling form [15J: 



where A is a transition rate of the original CP model, D is the diffusion rate, 




(1) 
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ranging between and 1, 6 is a critical exponent associated to the density variable 
g, corresponding to the value predicted by the mean-field approximation. The scaling 
function F(z) is singular at a point z = zq of its argument, such that the critical line, 
in the neighborhood of the multicritical point, corresponds to 



where <p is the crossover exponent, and the critical exponent of g on this line is 
determined by this singularity, being in general different from 9. 

One of the first studies of this problem was performed by Dickman and Jensen 
[H], who considered the model using supercritical series in A with the diffusion rate D 
taken as a fixed parameter. Therefore, in their calculation series expansions are derived 
for fixed values of the diffusion rate D, and analysing the series leads, among other 
information, to the phase diagram of the model with diffusion. However, they found that 
the fluctuations of the estimates provided by d-log Pade approximants show increasing 
fluctuations as the diffusion rate grows, so that the critical curve was obtained only up 
to D w 0.8. Since the crossover exponent characterizes the critical curve close to the 
infinite diffusion rate limit D —>■ oo no precise estimate of the crossover exponent was 
possible. The disapointing performance of the Pade approximants as the multicritical 
point is approached is not surprising, since it is known tha one-variable series analysis 
techniques are not effective close to such points [15]. More recently |10| . the model was 
simulated made in a the ensemble where the number of particles is conserved. These 
simulations display smaller fluctuations in the estimates, enabling the estimation of the 
critical line up to values near to the multicritiacal point, an leading the value <fi = 4.03(3) 
for the crossover exponent. 



This result is consistent with the lower bound > 1 predicted by Katori in [16] 



In the present work, we obtain estimates for the critical line and the crossover 
exponent <ft using a two-variable supercritical series analyzed using partial differential 
approximants (PDAs) [15J, [TFJ [18] . This technique seems to be more appropriate for the 
analysis of a two-variable series with a multicritical behavior, as shown by the results 
obtained from other models [Hi [19] . This paper is organized as follows. In section II we 
present the model and the mean-field results, in section III we show the derivation of 
the supercritical series and in section IV the analysis of this series is presented. Finally, 
in section V we the conclusions and final discussions of this work may be found. 

2. Definition of the model and mean-field results 

In a ci-dimensional lattice each site can be empty or occupied by a particle, so that 
we will associate an occupation variable rji — 0, 1 to the site i. The evolution of the 
system is governed by markovian local rules such that the particles are annihilated with 
rate 1 and created in empty sites with a transition rate Xn/z, where A is a positive 
parameter, n is the number of occupied nearest neighbors and z is the total number of 
nearest neighbors. In addition to these rules that define the CP, we include a diffusve 



(D c -D) = *b(A c -A)* 



(2) 
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Figure 1. Left: phase diagram obtained using the two-site mean-field approach. 
Right: the log-log plot of the same quantity is plotted and compared to the result 
<f> = 1 for this approach. 



process, allowing the hopping of particles to empty nearest neighbour sites at the rate 
D = D/(l — D). The configuration such that all sites are empty is an absorbing state. 
The transition from an active steady state, with a nonzero density of particles, to the 
absorbing state defines a transition line in the (A, D) plane as shown in figure [TJ 

A mean-field approach for this model can be obtained at several levels of 
approximation pp. In the one-site level the role of the diffusion is irrelevant since it 
contributed equally to creation and annihilation of particles at a given site i. Already 
in the two-site level it is possible to determine the critical line by using as variables the 
parameters A and D. This line is described by the expression 



A ' 2 



2- D+ V^ 2 + 4 



and it is not difficult to show that in the neighborhood of the multicritical point, 
(A c = 1, D = 1), the behavior of this curve is given by the scaling relation: 

(1 - £>) ~ (A - 1)* (4) 

where = 1. This crossover behavior may be seen in figure HJ 

This result is in accordance with the lower bound, <fi > 1, determined by Katori |16j . 
But, as is well known , the mean-field approach always overestimates the supercritical 
region of the models so that it is not surprising that more precise results for the exponent 
4> will differ from this unitary value. In the next section we will derive the supercritical 
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series in the variables A and D to determine the value of this exponent and compare it 
with that obtained in [TQl EH]- 



3. Derivation of the supercritical series for the model 

We use the operator formalism proposed by Dickman and Jensen [I] in order to derive a 
supercritical series. To this end, we define the microscopic configuration of the lattice, 
\r]), as the direct product of kets = <^i i \rji), with the following orthonormality 
property, (r]\r]') = Yii^Vi,Vi'- The particle creation and annihilation operators at site 
i are defined as 

Allrji) = (l-rj^rji + l), 

MVi) = VilVi - 1)- (5) 
In this formalism, the state of the system at time t is 

\m) = Y,p(v,t)\v), (6) 

{>?} 

where p(r], t) is the probability of a configuration r\ at time t. If we define the projection 
onto all possible states as ( | = S{r/}( ? ?l then the normalization of the state of the system 
may be expressed as ( l^} = 1- In this notation, the master equation for the evolution 
of the state is: 

The evolution operator S may be expressed in terms of the creation and annihilation 
operators as S = fiR + V where 

it! = DJ2(1 ~ 4-iA)A-i4 + (1 - A\ +1 Ai)A l+1 A\ + 

i 

+ J2(A i -AlA i ), (8) 

i 

V = £(4 - A,4)(4_ 1 A_ 1 + A\ +1 A i+1 ), (9) 

i 

where /i = 2/A. 

We notice that the operator R diffuses (01 — > 10) or annihilates particles (1 — ► 0), 
while the operator V acts in the opposite way, generating particles (0 — > 1). It is 
convenient to join the diffusion with the annihilation process to avoid ambiguities in 
the truncation of the series at a certain order. For small values of the parameter /i 
the creation of particles is favored, and the decomposition above is convenient for a 
supercritical perturbation expansion. Using the equations ([8]) and (Q the action of each 
operator on a generical configuration (C) is given by 



R{C) = D 



£(C0 + E(c; r ) + E(^ 



+ 



+ £(0 + -[(ri + 2r 2 ) + D + r](C), (i ) 



A two-variable series for the contact process with diffusion 



6 



where the first sum is over r\ sites with particles and one empty neighbor, the two next 
sums are over r 2 sites with particles and two empty neighbors and the last sum is over 
all sites occupied by a particle. Configuration (£■) is obtained moving the particle at 
the site i to the single empty neighbor site, {Cf r ^) is a configuration where the particle 
at the site i moved to the empty neighbor at the right (r) or at the left (/)is replaced 
by a hole and one of the empty neighbors (at the right or left). In the other hand, the 
action of operator V is 



V{C) = £(Cf ) + 2£(C;") - (gi + 2q 2 )(C), 



(11) 



where the first sum is over the q\ empty sites with one occupied neighbor, the second 
sum is over the q 2 empty sites with two occupied neighbors. Configuration (C"') is 
obtained occupying the site i in configuration (C). 

To obtain a supercritical expansion for the ultimate survival probability of particles, 
we start by remembering that in order to access the long-time behavior of a quantity, 
it is useful to consider its Laplace transform, 



m*)) = i e-«m))- 



(12) 



we find 



Inserting the formal solution \ip(t)) = e ^(O)) in the master equation 

Ms)) = (s-S)- 1 \m)- (13) 
The stationary state |^(oo)) = lim t ^oo {^(t)) ma Y then be found by means of the relation 

(14) 



l^(oo)) = lims|V>(s)}. 

A perturbative expansion may be obtained by assuming that 
powers of \i and using ( 1131 ) , 



(s)) can be expressed in 



$(s)> = l^0>+/i|^l>+/i 2 |^2> + 



V-fiR)- 1 ^)). 



Since 



we arrive at 



s-v-iiR)- 1 = (s-vy 1 h + fi( s - vy l R + fi 2 { s - vy 2 R 2 + 



l^o) 



vy l \m) 



(15) 



(16) 



;i7) 



and 



\^ n ) = (s-V)- 1 R\ij n . 1 ), (18) 

for n > 1. The action of the operator (s — V)~ l on an arbitrary configuration (C) may 
be found by noticing that 

{a - vy\C) = s - 1 {(c) + (s- vy l v{c)} , (19) 

and using the expression [TT] for the action of the operator V, we get 



vy\c) 



(C) + (s- vy 



E(C) + 2E(c 



7) 



(20) 
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where the first sum is over the q% empty sites and one occupied neighbor, the second sum 
is over the q 2 empty sites and two occupied neighbors, and we define s q = l/(s+qi + 2q 2 ). 

It is convenient to adopt as the initial configuration a translational invariant one 
with a single particle (periodic boundary conditions are used). Now, looking at the 
recursive expression ( 1201) . we may notice that the operator (s — V)^ 1 acting on any 
configuration generates an infinite set of configurations, and thus we are unable to 
calculate \ip) in a closed form. However, it is possible calculate the extinction probability 
p(s), which corresponds to the coefficient of the vacuum state |0). As happens also for 
models [H [3] related to the CP, configurations with more than j particles only contribute 
at orders higher than j, and since we are interested in the ultimate survival probability 
for particles P^ = 1 — lim s ^ sp(s), s q may be replaced by 1/q in equation fl20l) . An 
illustration of this procedure may be found in a previous calculation 

The algebraic operations described above is performed by a simple algorithm wich 
permit us to calculate 24 terms with a processing time to about 2 hours. Actually, 
the limiting factor in this operation is the memory required. In this way we define the 
coefficients b it j for the ultimate survival probability as: 



i i 24 i-2 

= 1 - -a - -a 2 - £ £ b,.yiV. (21) 

Z 4 i=3 j=0 



and they are listed in Table [Q 



4. Analysis of the series 

To obtain estimates of the critical properties, specially the critical line, from the 
supercritical series for the ultimate survival probability as given by the equation (l2Tj) . 
we initially use d-log Pade approximants. These approximants are defined as ratios of 
two polynomials 



Qm(\) i • Ejf , <ij\-> 



Flm(X) = = ; r< = fW- ( 22 ) 



In our case the function /(A) represents the series for 4- lnPoo(A). As /(A) is a function 
of one variable, we fix the value of D to calculate these approximants. For a fixed 
value of D one pole of the approximant F will correspond to the critical point while 
the associated residue will be the critical exponent (3. We calculate approximants with 
L + M < 24, restricting our calculation to diagonal or close to diagonal approximants, 
which usually display a better convergence. Thus L = M + £, where £ = 0, ±1 and 
with D — D/(l + D) ranging between and 0.8. Examples of estimates for the critical 
values of fi obtained from these approximants is given in Table [2] for different values of 
the diffusion. 

For each value of the diffusion rate, we calculate about eight approximants, 
obtaining the estimate of /i c associated to diffusion as an arithmetic average of results 
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Table 1. Coefficients for the series expansion for ultimate survival probability 
corresponding to the CP model with diffusion. The indexes refer them to the equation 

(ED- 



'hi 



'■J 



4 



6 



7 




1 


1 

2 

1 
2 
3 

1 
2 
3 
4 

1 
2 
3 
4 
5 

1 
2 
3 
4 
5 
6 

1 
2 
3 
4 
5 
6 
7 

10 
1 

2 
3 
4 
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0.25000000000000000000 x 10° 
-0.25000000000000000000 x 10° 
0.28125000000000000000 x 10° 
-0.37500000000000000000 x 10° 
0.37500000000000000000 x 10° 
0.34375000000000000000 x 10° 
-0.50781250000000000000x10° 
0.57812500000000000000x10° 
-0.62500000000000000000 x 10° 
0.44726562500000000000 x 10° 
-0.76220703125000000000x10° 
0.85058593750000000000 x 10° 
-0.83984375000000000000 x 10° 
0.10937500000000000000x10° 
0.60223388671874955591 x 10° 
-0.11734619140625004441 xlO 1 
0. 15190429687499997780 x 10 1 
-0.14140624999999984457X10 1 
0.10878906249999982236X10 1 
-0.19687500000000000000X10 1 
0.83485031127929687500x10° 
-0. 18110389709472716202 x 10 1 
0.25234603881835981909 x 10 1 
-0.29291381835937464473 x 10 1 
0.24864501953125062172X10 1 
-0.10754394531250017764x 10 1 
0.36093750000000000000 x 10 1 
0.11814667913648828623X10 1 
-0.28569926950666579835 x 10 1 
0.42781094621729156557x 10 1 
-0.48761836864330092567x 10 1 
0.54410674483687788694 xlO 1 
-0.48496839735243062464 xlO 1 
0.11848958333333414750x10° 
-0.67031250000000000000 xlO 1 
0.16988672076919952847X10 1 
-0.45030223008843064392 x 10 1 
0.73700355965291270977x 10 1 
-0.92486491500105252328 x 10 1 
0.87502182305104447835 x 10 1 



11 



12 



13 



6 
7 
8 

1 
2 
3 
4 
5 
6 
7 



14 




1 

2 
3 
4 
5 
6 
7 
8 
9 

10 



1 

2 
3 
4 
5 
6 
7 
8 
9 

10 
11 

1 

2 



0.10567643059624561630xl0 2 
0.34998474121093847700x 10 1 
0.12568359375000000000x 10 2 
0.24775957118666096513x 10 1 
-0.71703082845177412707x 10 1 
0.12429764028158224676x 10 2 
-0.16676337106809967281 xlO 2 
0.19435605471721252968x 10 2 
-0.15230894658300556443x 10 2 
0.13923689787279895924x 10 2 
-0.24910518081099844778x 10 2 
-0.13853664539478481643x 10 2 
-0.23740234375000000000 x 10 2 

0.36488812342264926869 x 10 1 
-0.11443929729648042226x 10 2 
0.21418735868689896762x 10 2 
-0.29831307350977681381 x 10 2 
0.34272964785342651339x 10 2 
-0.40398672142671166796x 10 2 
0.27855684964608855125x 10 2 
-0.13595902518316915319x 10 2 
0.60935236387946176251 x 10 2 
0.40270442479922508028 x 10 2 
0.45106445312500000000x 10 2 
0.54293656084851154020x 10 1 
-0.18322144692814863021 x 10 2 
0.36259195896082665911 xlO 2 
-0.54866931326313050477x 10 2 
0.67130818799164941879x 10 2 
-0.64293858561553619779x 10 2 
0.81638245065997452343x 10 2 
-0.59206203158103356543x 10 2 
-0.11386872839957611347X 10 2 
-0.15017302299755911577x 10 3 
-0.10358591484729181786x 10 3 
-0.8611230468750000000x 10 2 
0.8132542219307161701600x 10 1 
-0.29467694610727896531 x 10 2 
0.62075441392908530247x 10 2 
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j J hi 

14 4 0.12740413805065173847xl0 3 

5 -0.14872799806680012580xl0 3 

6 0.11055259565281477308xl0 3 

7 -0.15322889380857196784xl0 3 

8 0. 15125687744198603468 xlO 3 

9 0. 12100775973170790678 xlO 3 

10 0.36717258144235222517xl0 3 

11 0.24952145042880511028xl0 3 

12 0.16504858398437500000 xlO 3 

15 0.12275012836144505002 xlO 2 

1 -0.47363165128788978109xl0 2 

2 0.10546586796137503939 xlO 3 

3 -0.1762739398241103288xl0 3 

4 0.23298609118631188153xl0 3 

5 -0.27071715385838172097xl0 3 

6 0.33267266081610591755xl0 3 

7 -0.18014432368848466126xl0 3 

8 0.24323790718884771422 xlO 3 

9 -0.43351480900655008099xl0 3 

10 -0.48110863117407200207xl0 3 

11 -0.88513499744041246231 xlO 3 

12 -0.57707323452079549497xl0 3 

13 -0.31740112304687500000xl0 3 

16 0.18620961415130427241xl0 2 

1 -0.76547748518027589171 xlO 2 

2 0.17936794520034777634 xlO 3 

3 -0.30967915791812674797xl0 3 

4 0.45127497217248452444xl0 3 

5 -0.53698491310493750461 xlO 3 

6 0.51912309945306844838 xlO 3 

7 -0.74776029427868388666 xlO 3 

8 0.31390176074091152714xl0 3 

9 -0.23066668763466492464xl0 3 

10 0.12806582574012411442xl0 4 

11 0.15251186818533849419xl0 4 

12 0.21006816788719602300xl0 4 

13 0.12983815231244370807xl0 4 

14 0.61213073730468750000 xlO 3 

17 0.28405733950686048672xl0 2 

1 -0.12342415559750365617xl0 3 

2 0.30541526863109334045 xlO 3 



1 J 

3 -0.55737662171810006839X10 3 

4 0.80584299514354984240xl0 3 

5 -0.10663609887965685630 xlO 4 

6 0.12359653482806177180xl0 4 

7 -0.86967640846259655518 xlO 3 

8 0.16397080915531578285xl0 4 

9 -0.72176857313912660175xl0 3 

10 -0.39467789553376610456 xlO 3 

11 -0.37241724510229601037xl0 4 

12 -0.43433579993156563432 xlO 4 

13 -0.49111545112523144780xl0 4 

14 -0.28643360597728060384 xlO 4 

15 -0.11834527587890625000xl0 4 

18 0.435207828742268674200 xlO 2 

1 -0.20009555228747112210xl0 3 

2 0.51666085550451919062xl0 3 

3 -0.98221829678260564833 xlO 3 

4 0.15468813789798582548xl0 4 

5 -0.19452358996876919264xl0 4 

6 0.23142581042234874076xl0 4 

7 -0.29269143532222028625 xlO 4 

8 0.11929249112512670763xl0 4 

9 -0.33469260525748682085 xlO 4 

10 0.22662397929431922421 xlO 4 

11 0.33534053058169365613xl0 4 

12 0.10497352179518215053xl0 5 

13 0.11603782689493993530xl0 5 

14 0.11326231062527707763xl0 5 

15 0.62269615466431168898xl0 4 

16 0.22929397201538085938xl0 4 

19 0.66930218067969633466xl0 2 

1 -0.32354897975245813768 xlO 3 

2 0.88042629554806353553xl0 3 

3 -0.17413898514485581472xl0 4 

4 0.275 74673463423659996 xlO 4 

5 -0.39760107863198809355xl0 4 

6 0.45478180647223589403 xlO 4 

7 -0.44829538123020120111xl0 4 

8 0.71439950475520599866xl0 4 

9 -0.11433220550468702186xl0 4 

10 0.58699584417021997069xl0 4 
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% 


j 


hi 


i 


3 


hi 


19 


11 


-0.80007309950477119855xl0 4 




15 


-0.19550932821 134978440 xlO 6 




12 


-0.13849586891129882133x 10 5 




16 


-0.17865808405461237999xl0 6 




13 


-0.28630836746692133602x 10 5 




17 


-0.13001245697926016874xl0 6 




14 


-0.29699369329023520550x 10 5 




18 


-0.60334798749708410469 x 10 5 




15 


-0.25808082841629722679 x 10 5 




19 


-0.16853935146331787109x 10 5 




16 


-0.13385541065918856475x 10 5 


22 





0.24876519640902955643xl0 3 




17 


-0.44510006332397460938x 10 4 




1 


-0. 13849806980532957823 x 10 4 


20 





0.10337399908883011790x 10 3 




2 


0.42491525612070690840X 10 4 




1 


-0.52548922262251915072x 10 3 




3 


-0.95509490156905540061 x 10 4 




2 


0.14837268484319015442x 10 4 




4 


0.17072122607398174296xl0 5 




3 


-0.30794482965317188246xl0 4 




5 


-0.25259991312968326383x 10 5 




4 


0.51825134723219762236x 10 4 




6 


0.36959825267281092238xl0 5 




5 


-0.70384028435684167562 x 10 4 




7 


-0.39724946702482979163 x 10 5 




6 


0.95392163017937873519x 10 4 




8 


0.4108248820168149541 1 x 10 5 




7 


-0.10775090683332531626x 10 5 




9 


-0.69996295461125846487x 10 5 




8 


0.72306238600702890835x 10 4 




10 


-0.47433844045513687888X 10 4 




9 


-0.17556031650766155508x 10 5 




11 


-0.94472082453477501986 x 10 5 




10 


0.45576565336850308086 x 10 3 




12 


0.18231631734823597071 xlO 5 




11 


-0.68853229355527937514x 10 4 




13 


0.74507403911089320900 x 10 5 




12 


0.27725591978402942914x 10 5 




14 


0.28127919612702319864x 10 6 




13 


0.46620453043310422800x 10 5 




15 


0.40344708716105029453x 10 6 




14 


0. 75787894589888033806 x 10 5 




16 


0.49346542769156675786x 10 6 




15 


0.73699439769879478263x 10 5 




17 


0.42523410507562680868x 10 6 




16 


0.58191226154708187096xl0 5 




18 


0.28817317888963740552x 10 6 




17 


0.28519783989193914749 x 10 5 




19 


0. 12690253004534545471 xlO 6 




18 


0.86547234535217285156 xlO 4 




20 


0.32865173535346984863X 10 5 


21 





0. 15998830271612598608 x 10 3 


23 





0.38696067609131841891 x 10 3 




1 


-0.85206464313625656359x 10 3 




1 


-0.22521288931067451813 xlO 4 




2 


0.25270411365871177622xl0 4 




2 


0.72239403631839231821 x 10 4 




3 


-0.53883333002320678133X 10 4 




3 


-0. 16540644823073282168 x 10 5 




4 


0.93069165541648508224x 10 4 




4 


0.30965978110985601234xl0 5 




5 


-0. 14256926805592749588 x 10 5 




5 


-0.49715203648917136888xl0 5 




6 


0. 16894039347875652311 x 10 5 




6 


0.63140021723358484451 x 10 5 




7 


-0.21057424914324066776 x 10 5 




7 


-0.90160716109902248718xl0 5 




8 


0.2676311383571557235xl0 5 




8 


0.96095356087432839558 x 10 5 




9 


-0.76173492536429366737x 10 4 




9 


-0.63690648044614848914xl0 5 




10 


0.42006937485571783327x 10 5 




10 


0.18830794813510467065x 10 6 




11 


-0.12698208721358355433x 10 4 




11 


0.57163083968398816069x 10 5 




12 


-0.53474107987358220271 x 10 4 




12 


0.19078924010583921336x 10 6 




13 


-0.90835383726464555366 x 10 5 




13 


-0.11462564856718564988xl0 6 




14 


-0.14158471202173284837x 10 6 




14 


-0.34295621170633088332xl0 6 
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i 


3 




23 


15 


-0.83015840201854216866x 10 6 




16 


-0.11002212525564364623x 10 7 




17 


-0.12227985275158903096x 10 7 




18 


-0.99733247054306510836 x 10 6 




19 


-0.63431234698974736966 x 10 6 




20 


-0.26563714369463571347x 10 6 




21 


-0.64165338807344436646 x 10 5 


24 





0.60509199550873199769x 10 3 




1 


-0.36606232348859748527x 10 4 




2 


0.12139375556095224965xl0 5 




3 


-0.29330923818458009919x 10 5 




4 


0.55510496259245075635 x 10 5 




5 


-0.89523794399456470273x 10 5 




6 


0.13690056986421268084x 10 6 




7 

8 


-0.14950913041347730905x 10 6 
0. 20276027262469305424 xlO 6 




9 


-0.24868958310934680048x 10 6 




10 


0.44658906685524416389 x 10 5 




11 


-0.50473239786133670714x 10 6 




12 


-0.20583143731838543317x 10 6 




13 


-0.31319042987920023734x 10 6 




14 


0.51903724071582528995 x 10 6 




15 


0.12406518345377091318xl0 7 




16 


0.23558257104359627701 x 10 7 




1 ( 


0.29049027840298512019x10 




18 


0.29832956761808455922 x 10 7 




19 


0.23110049545479607768 x 10 7 




20 


0. 13877345160856433213 x 10 7 




21 


0.55381177327087428421 x 10 6 




22 


0.12541407130479812622xl0 6 



funnished by these set of approximants. From this we obtain the phase diagram shown 
in the figure [2j With the purpose of comparison with the results coming from the 
conservative simulations [THj we use the variables a = fi/2 and -D e // — aD/{l + aD). 

For higher values of the diffusion the dispersion increases, and estimates with larger 
error bars are found. Nevertheless, the exponent = 4 seems to describe well the 
calculated points of the critical line. However, this result would be different if we 
used approximants to series closer to the infinite diffusion rate limit. This error in 
the approximants for high values of the diffusion rate is attributed to the alternated 
sign of the series terms [H] . Another explanation would come from the fact that in the 
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\j 


T 

ij 


M 

1V1 


fJ-c 


n 
U 


1 n 
1U 


1 n 

1U 


U.0U040 




1 1 
11 


1 1 


U.DU04D 


n i 


i n 


i n 






n 


11 


0.62266 


0.7 


10 


10 


0.85353 




11 


11 


0.84513 


0.8 


10 


10 


0.96256 




11 


11 


0.94246 



Table 2. Esimates for ritical points obtained by d-log Pade approximants. Note that 
as the value of I? = D/(l + D) grows the dispersion in the value estimates also grows. 




Figure 2. At left panel we have the phase diagram obtained by two-site mean-site 
approach. At right the log-log plot of the same quantity is plotted showing that value 
of the crossover exponent (f> is unity in this approach. 



neighborhood of a multicritical point the reduction of a two- variable series to one variable 
leads to very poor estimates of the critical properties pj] close to a multicritical point. 
To overcome this problem we analyze the series using Partial Differential Approximants 
(PDA's) [15], that generalize the d-log Pade approximants for a two- variable series. 
These approximants are defined by the following equation 

P L (x, y)F(x, y) = Q M (x, y)^^ + J2n(*, v)^^, (23) 
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where P, Q, and R are polynomials in the variables x and y with the set of nonzero 
coefficients L, M, and N, respectively. The coefficients of the polynomials are obtained 
by substitution of the series expansion of the quantity which is going to be analyzed 

f(x,y)= £ f(k,k')x k y k ' (24) 

k,k'=0 

into equation ([23]) and requiring the equality to hold for a set of indexes defined as 
K. This procedure leads to a system of linear equations for the coefficients of the 
polynomials, and since the coefficients f^k 1 of the series are known for a finite set of 
indexes this places an upper limit to the number of coefficients in the polynomials. 
Since the number of equations has to match the number of unknown coefficients, the 
numbers of elements in each set must satisfy K = L + M + N — 1 (one coefficient is 
fixed arbitrarily). An additional issue, which is not present in the one- variable case, is 
the symmetry of the polynomials. Two frequently used options are the triangular and 
the rectangular arrays of coefficients. The choice of these symmetries may be related to 
the symmetry of the series itself [PTj . 

It is possible to show [17] that we can determine the multicritical properties using 
the equation [23] and the hypothesis of that in the neighborhood of the multicritical 
point, the function / obeys the following scaling form 

f(x,y)*\AxrZ (25) 



where 



and 



|Ax| 



Ax = (x - x c ) - (y - y c )/e 2 , (26) 



Ay = (y - y c ) - ei(x - x c ). (27) 

Here v is the critical exponent of the quantity described by / when Ay — 0, e\ and e 2 
are the scaling slopes [TJ] and is the crossover exponent. 

On the other side, our calculation was successful when we use the method of the 
characteristics to integrate equation ( [23]) . This is made by introducing a timelike variable 
r, so that a family of curves is obtained in the plane (rr(r), y(r)) (the characteristics). 
Such curves obey to the equations 

-j- = QM{x{r),y{r)), 

^- = R N (x(r),y(r)). (28) 
dr 

It is possible to show that integrating the equations ([281) from a specific point of the 
critical line, the resulting characteristic is equivalent to the the critical line. In figure 
E] we show a comparision between a characteristic curve and the simulational result 
[TO] . The number of elements in the sets of the calculated approximants was varying as 
follows: K = 55 - 190, M = N = 20 - 53 and L = 15 - 54. 
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In each of these curves, we calculate his inclination in the neighborhood of the 
multicritical point, determining a value for the exponent <fi and the mean value for this 
exponent results as <fi = 4.02±0.13, consistent with simulations in the particle conserving 
ensemble |10j . 




Figure 3. Comparision between a characteristic curve (solid line) and a numerical 
simulation result (circles). The coincidence is evident, including the region of the 
infinite diffusion limit (inset). 

Using all the characteristic curves calculated we derive an 'average curve', 
calculating for each value of a on arithmetic average for D e ff. This curve is shown 
in the figure H] jointly with the result originating from the simulation and with the 
scaling form (1 — D e ff) ~ (1 — a) 4 . In the same figure we see that the exponent <fi = 4 is 
well fitted to the results of simulation and of the PDA's. This scaling form is based on 
the argument of the scaling function Z presented in the equation (I25j) . where = 4 and 
zq is a parameter properly chosen. We remark that this scaling form coincides with the 
characteristic curve and with the simulation even in the weak diffusion regime. This is 
somewhat surprising since its validity would be expected only in the neighborhood of 
the multicritical point. 

Unfortunately, even using the algorithm proposed by Styer [17] we were not able 
to obtain precise estimates for the crossover exponent </> from the scaling form shown 
in equation (1251) . However, integrating a set of approximants, we could determine the 
characteristic curves whose initial point is conicident with the critical point of the CP 
without diffusion of particles. These curves are estimates for the critical line of the CP 
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-2 -1.5 -1 -0.5 

In (1-a) 



Figure 4. Comparision between a characteristic 'average curve' (solid line), the 
numerical simulation result (triangles) and the curve obtained from the scaling form 
(1 — D e ff) ~ (1 — a) 4 (circles). At the right panel we see that phi = 4 looks a good 
estimative for the crossover exponent. 

model with diffusion going beyond the values achieved in [13] and [10] and corroborating 
the initial result of this last reference in that 0^4. 

5. Conclusion 

Calculating a supercritical series for the ultimate survival probability and analysing it 
using PDA's we obtain estimates for the critical line of the CP model with diffusion. The 
critical line was derived through the integration of the equation ( 1231) by the method of 
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the characteristics. Direct results for the value of the crossover exponent using the 
scaling form [25] using Styer's algorithm [IT] were not possible to be obtained with 
an acceptable precision. However the method of the characteristics permitted us to 
calculate the critical line and, consequently, the value for the crossover exponent 0. Our 
result, <p = 4.02 ± 0.13 is in agreement with that derived in pU] and explores a region 
of diffusion very close to the multicritical point. 

The technique of the two- variable supercritical series associated with PDA analysis 
was shown to be accurate enough to determine the critical properties in similar models 
[Tlj [20] . Therefore, we believe that a natural extension for this work is analyze related 
models that apparently possess non-trivial multicritical points. This semms to be the 
cases of the pair-creation and triplet-creation models with diffusion, also studied in [10J. 
This research is already in course. 

W.G. Dantas acknowledges the financial support from Fundagao de Amparo a 
Pesquisa do Estado de Sao Paulo (FAPESP) under Grant No. 05/04459-1. 
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